Large anisotropy in the optical conductivity of YNi2B2C 
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The optical properties of YNi2B2C are studied by using the 
first-principles full-potential linearized augmented plane wave 
(FLAPW) method within the local density approximation. 
Anisotropic behavior is obtained in the optical conductivity, 
even though the electronic structure shows 3D character. A 
large peak in is obtained at 2.4 eV. The anisotropic opti- 
cal properties are analyzed in terms of interband transitions 
between energy levels and found that the Ni site plays an im- 
portant role. The electronic energy loss spectroscopy (EELS) 
spectra are also calculated to help elucidate the anisotropic 
properties in this system. 
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Quaternary transition metal borocarbide supercon- 
ductors with formula RM2B2C (R= Y, rare-earths; 
M=transition metal) have attracted much attention due 
to their relatively high superconducting transition tem- 
peratures in the intermetallic compounds including ex- 
change enhanced transition metal elements such as Ni, 
Pd, and Ru QJ^. Of special interest is the coexis- 
tence and interplay between superconductivity and mag- 
netism in the systems containing rare-earth elements . 
YNi2B2C with Tc—\b.Q K serves as a reference system 
for understanding the superconducting mechanism in this 
family, because there are no complications introduced by 
the presence of magnetic rare-earth elements. 

YNi2B2C crystallizes in a body-centered tetragonal 
structure with a ThCr2Si2-type (space group lA/mmm, 
Z?4^). The crystal structure is layer- like in the c axis, 
reminiscent of the high- Tc cuprate superconductors, con- 
sisting of two alternating layers of Ni2B2 and YC. Elec- 
tronic structure studies show three dimensional (3D) 
character and so it is thought to be a conventional BCS 
type superconductor with a relatively high density of 
states (DOS) at the Fermi level 1^-0]. Despite the layered 
anisotropic crystal structure, some physical properties 
of YNi2B2C are often considered to be isotropic. Note 
that the 3D character of the electronic structure does 
not necessarily mean the isotropy of physical properties. 
In fact, it has not yet been settled whether the physi- 
cal properties of YNi2B2C are isotropic or not. Civale 
et al. 1^ observed an anisotropic effective mass ratio 
7 ~ ("ima2;/"imm)^^^ l-l; whilc Johnston-Halpcriu 
et al. PI observed an isotropic result, 7 ^ 1.005. Fisher 
et al. |lC|] reported isotropic resistivity for non-magnetic 
YNi2B2C and anisotropic resistivity for magnetic boro- 



carbide superconductors like ErNi2B2C and HoNi2B2C. 
Rathnayaka et al. [ p| observed an isotropic upper criti- 
cal field, Hc2, for YNi2B2C, but a small anisotropic Hc2 
for LaNi2B2C. In contrast, anisotropic properties were 
reported for YNi2B2C thin-films in the upper critical 
field , thermal conductivity |]l3| , critical field for vor- 
tex lattice transition Jl^ , and paramagnetic susceptibil- 
ity from nuclear magnetic resonance (NMR) experiments 
JlSf . von Lips et al. ||l^] used X-ray absorption spec- 
troscopy (XAS) and the linear combination of atomic-like 
orbitals (LCAO) band method for YNi2B2C to determine 
the anisotropic nature of the unoccupied DOS. 




FIG. 1. (a) Crystal structure and (b) the first Brillouin 
zone of YNiaBaC. 

The optical properties of YNi2B2C have been studied 
several times since the first measurement by Widder et al. 

, who obtained anisotropic results by linking together 
optical reflectance and electron energy loss spectroscopy 
(EELS) measurements. On the other hand, Bommeli et 
al. ||l^ reported that the structural anisotropy would not 
affect the optical data based on the observation that both 
the polycrystalline and single crystal samples produced 
equivalent experimental results. Therefore, the results 
were conflicting, and until recently, the anisotropy prob- 
lem in the optical properties of Ni borocarbides has not 
been resolved [^^-^ . It is important to understand cor- 
rect optical properties since the superconducting energy 
gap is obtained from optical measurements pq] . 
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Interestingly, all the above optical groups utilize the 
3D electronic structures of Ni borocarbides for their ex- 
periments. It thus indicates that the electronic structure 
information only may not be enough to determine def- 
initely the optical anisotropy. In this Letter, to resolve 
this controversial issue, the optical conductivity and the 
electron energy loss spectroscopy (EELS) have been stud- 
ied theoretically for the first time by using first principles 
calculations. We have found that the optical properties 
of YNi2B2C are really anisotropic. 

Calculations are performed by employing the highly 
precise first principles local density full-potential lin- 
earized augmented plane wave (FLAPW) method [^ , 
which makes no shape approximations to the charge den- 
sity or the potential. The Hedin-Lundqvist form |^3| is 
employed for the LDA exchange-correlation energy. We 
used lattice parameters determined by experiment [p4[ , 
240 special k-points ||2^ for the self-consistency cycle and 
8619 k-points for the calculation of the optical proper- 
ties inside the irreducible Brillouin zone. The optical 
conductivity arises from the optical transitions of the oc- 
cupied valence electrons to unoccupied states above Ej? 
pGf . Each interband transition is broadened in the cal- 
culation by introducing a finite phenomenological inverse 
life time, A=0.05 eV. The Drude term due to free carriers 
is taken into account by using parameters from experi- 
ment 0. 

Figure |l| shows the crystal structure and the first Bril- 
louin zone of YNi2B2C. In Fig. |l](b), points a,b,c, and 
d are not symmetry points since the symmetry is not a 
local maximum at those points. Points a and d (and b 
and c likewise) are the same points by the translational 
symmetry in k-space. X has the highest symmetry as F 
with 16 symmetry elements. 

Figure ^ shows the calculated energy band structures 
along the symmetry lines. Band dispersions agree well 
with existing calculations. However, in the present full- 
potential band calculation, the 19th band at F, which 
gives a small calabash shaped Fermi surface in the 
LMTO-ASA calculation is lifted up to higher energy. 

Figure ^ compares the calculated optical conductiv- 
ities for YNi2B2C with experiments, in which and 
az represent theoretical conductivities along the x and 
z directions, respectively. One can see a clear difference 
between and az, implying the anisotropic nature of 
the optical conductivity. Anisotropic optical conductivi- 
ties in the calculation are consistent with the anisotropy 
observed in the experiment by Widder et al. [ p^ , repre- 
sented by CTj^-w and cr^.w in Fig. ^. Peaks are located 
at 2.48, 7.6, 8.9, and 9.6 eV for cr^, and at 2.4, 5.0, 8.16, 
and 9.5 eV for tr^ in the calculation. The most prominent 
and anisotropic peak in the calculation is the one at 2.4 
eV of az, while a.^ has a peak at 2.48 eV with a smaller 
magnitude than cr^. However, the intensity and position 
of the peak of az are somewhat different from those in 
experiment: the highest peak is located at 1.6 eV, which 



is lower by 0.8 eV than the theoretical one. It is unusual 
that the experimental peak is located at lower frequency 
considering that, in most LDA calculations, theoretical 
peaks are found at lower frequencies than experimental 
ones. The disagreement is to be clarified further theoret- 
ically or experimentally. Other optical properties such 
as the reflectivity R and the dielectric constants e show 
similar anisotropics which will be discussed elsewhere. 




ZbN Pa Zr X 



FIG. 2. Electronic energy band structure of YNi2B2C. 
Dot and solid arrows indicate typical interband transitions 
along the x and z directions, respectively (see text and Fig. 
0). 

Bommeli et al. obtained isotropic optical conduc- 
tivity (denoted by ax-h in Fig. ||). The overall shape 
of a^-h is seemingly close to a^-w, except that some 
fluctuations occur when the relative magnitude of (Tz_w 
with respect to ax-w becomes large. This feature can 
be understood if a^-h corresponds to the average of the 
anisotropic optical conductivity. Although it is not clear 
in ax-h, there is a peak at around the 2.4 eV in ax-w. 
The 2.4 eV peak in ax is also observed in a recent exper- 
iment by Lee et al. pO| ] . Hence the experimental peak at 
2.4 eV can be assigned to the calculated ax peak at 2.48 
eV, 

The optical conductivity is obtained directly from the 
optical transitions between energy levels in the solid. 
Hence the peaks in the optical conductivity can be as- 
signed to specific interband transitions. However, peaks 
in derived quantities from the optical conductivity, such 
as the reflectivity R and the refractive index n, may be 
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shifted from the peaks of the optical conductivity. In 
order to identify the peaks manifested in the optical con- 
ductivity, we plot in Fig. ||the momentum resolved opti- 
cal conductivity cr(k,uj) along the symmetry lines in the 
Brillouin zone. Integrating (T(k, lu) over k-space gives rise 
to the optical conductivity a{uj). 

Figures ^ (a) and (b) provide the optical transitions 
along the x and z directions, respectively. The darker in- 
tensity in Fig. |4| represents the stronger interband tran- 
sitions. It is evident from Figs, ^(a) and (b) that the 
optical transitions are anisotropic between the x and z 
directions. The anisotropy is due to different matrix ele- 
ment between the x and z directions. Strong transitions 
in the z-direction occur around P and N, while they oc- 
cur at points near F in the x-direction. The anisotropy 
is most prominent at the 2.4 eV transition near P; there 
is a strong transition in the z-direction, but a negligible 
transition in the x-direction. The intensity of the 2.4 
eV transition is further enhanced by the transition near 
N. In Fig. H, corresponding transitions for 2.4 eV are 
marked as small solid arrows at N and P. Other strong 
transitions identified in Fig. ^ are also marked as arrows 
in Fig. 1^. It is noteworthy that the unoccupied bands 
(19th band) for the 2.4 eV transition near P have similar 
shapes to the occupied bands (18th band), reminiscent 
of the nesting of bands. 
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FIG. 3. Optical conductivity of YNi2B2C. (Jx and rep- 
resent results from the calculation. Experimental data are 
taken from Ref. [17](aa:-w and a^-w) and Ref. [18](cri:-b). 

In order to see the detailed nature of the transition, we 



have calculated contributions from each atomic species 
for the bands related to the 2.4 eV transition at P. The 
18th band consists of Ni-d322_r2 (36%) and Ni-dxy{55%) 
states, while the unoccupied band (19th band) consists of 
Ni-d^y(16%), Ni-s(12%), Ni-]5^(4.4%), B-p(6.9%), and C- 
p(16%) states. Since the occupied band originates mainly 
from Ni-d states, the transition matrix element is deter- 
mined by Ni site. By the selection rule about the sym- 
metry of the wave functions, one can assign the 2.4 eV 
peak in az to the transition between Ni-d322_^2 and Ni- 
Pz- The transition matrix element for is zero in this 
case, which explains the negligible intensity at 2.4 eV 
near P in (T2.(k,a;). Thus Ni states play a vital role in 
yielding a large anisotropic peak in cr^. This is indeed 
consistent with the dominant Ni-d DOS near the Fermi 
level §. 




FIG. 4. Momentum resolved optical conductivity (j(k,u}) 
for (a) X and (b) z direction, respectively. 

Figure ^ depicts the loss function, Im(— l/e(ti;)), of 
the EELS measurement with zero momentum transfer. 
There is no EELS data available for single crystals. 
Hence, in the figure, the calculated result is compared 
with Widder et a/.'s experiment for polycrystals [ pT[ . 
There is no plasma peak due to free carriers which is ex- 
pected to occur between 3.75 eV and 4.25 eV 0,|l|,|§. 
It is damped due to interband transitions (Landau damp- 
ing), in agreement with the experimental EELS result. 
The calculated loss function shows a two peak structure 
in both the x and z-direction, and the spectrum in the 
z-direction is narrower and higher than that of x direc- 
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tion. Two peaks are located at 22.7 cV and 26 cV in the 
z-direction. In YNi2B2C, there are 33 valence electrons 
per formula unit. The simple plasma frequency relation, 
LUp = 47rne^/?7i, gives 26.3 eV which agrees well with the 
above plasma frequency of 26 eV. Therefore, one can ex- 
plain the plasma frequency at 26 eV as originating from 
the response of all valence electrons to the external elec- 
tric field. Then the peak at 22.7 eV can be ascribed to 
the valence electrons near Ep: Ni-rf, Y-d, B-p, and C-p 
states. In the experiment, two peaks arc observed at ~ 
20 eV and ~ 35 eV. Although the two peak structure 
in the calculation agrees qualitatively with experiment, 
the peak positions are different. The EELS experiment 
with single crystal is required to clarify the discrepancy 
between theory and experiment. 

In conclusion, wc have investigated the optical prop- 
erties of YNi2B2C by using the first-principles FLAPW 
method. The optical conductivity is found to be 
anisotropic between the x and the z direction. A strong 
peak at 2.4 eV is obtained in the z direction, which is as- 
signed to the transition near P between Ni-dsz^-r^ and 
Ni-p^ states, indicating that the Ni-site plays an impor- 
tant role in the anisotropy of YNi2B2C. The calculated 
EELS spectra are also found to be anisotropic with a two 
peak structure. 
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FIG. 5. Electron energy loss (EELS) spectrum Im(— 1/e) 
of YNi2B2C. Crosses represent experimental data reproduced 
from the experimental Widder et al. [17]. 
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